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We study the properties of a quantum impurity embedded in a superconducting host. The superconductor 
is described by the negative U Hubbard model while the impurity introduces a repulsive interaction to the 
system. We discuss the influence of this 'on-impurity' Coulomb repulsion on the local properties (density 
of states, electron pairing) of the superconductor. We show the condition of n - like behaviour, defined as 
two subsystems having a phase difference of n, in the system by using a proper combination of attractive 
pairing interaction and repulsive one located at impurity site. 
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1 Introduction 

The role of impurities in superconducting hosts has been a subject of intensive theoretical and experimental 
studies El El El El El recently. Mainly, the bulk properties, like transition temperature or density of states 
have been investigated 1 6 1 . However in a conventional (as discussed in this paper) s-wave superconductors 
(SC), non-magnetic impurities have only little effect on the transition temperature according the Anderson 
theorem ,7]|8l|9). On the other hand very little attention has been payed to the understanding the properties 
of the local environment of a single impurity, like local density of states (DOS) or spatial variation of 
the order parameter (see however 1 10 1). The effect of non-magnetic impurities on the local properties of 
superconductors with exotic pairing has been studied more extensively in 1 1 1 1 fl 2 1 and references therein. 
The large interest in the subject of impurities in exotic SC comes from the fact that the Anderson theorem 
does not work in the later case and impurities have more drastic influence on the properties of the high-T c 
and other non-s-wave superconductors. In particular, it is widely accepted that high-T c materials have 
d-wave pairing state 1 13 1, which means that SC order parameter changes its sign under ir/2 rotation. This 
leads to the so called 7r-phase behaviour and can be seen in nonuniform systems, like those with surfaces, 
vortices, cracks, twin boundaries or impurities 1 141 1 1 51 . 

A system is regarded as in the ir -phase if there is a sign change of the order parameter between two 
subsystems. The simplest example is the junction made from two superconductors with the phase of the 
order parameter equal to tt 1141 1 1 51 . In this case Josephson current becomes negative in contrast to the 
usual 0-phase junction. An other example are granular high-T c materials which can likely form network of 
microscopic ir -junctions 1 16 1 between small regions with different phases of the order parameter. In such 
systems the zero-energy Andreev bound states, zero-bias conductance peaks, paramagnetic Meisner effect 
and spontaneously generated currents take place 1 1411 1 51 . 

The situation is very similar to those in which paramagnetic impurities with strong on-site Coulomb 
repulsion are placed in classical s-wave superconductor. Due to the proximity effect |17jj the SC order 
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parameter is created on the impurity site. Moreover, because of the fact that the impurity coupling constant 
has opposite sign to SC one, there is a sign change of the order parameter between impurity and surrounded 
SC. Therefore, we regard the system as in 7r-phase. It is the purpose of the present work to clarify if such 
system has all necessary ingredients to regard its as a 7r-phase state. 



2 The model 

The system is described by negative U Hubbard model 1 18 1 with the Hamiltonian: 



H = ^2 (*« _ M<5y) 4a c i<y + 2 XI Ui n i<rni-<7, (1) 

ijcr icr 

where i, j label sites of a square lattice, fcy = — t is the hopping integral between nearest neighbour 
sites and (i is the chemical potential. Ui < describes attraction between electrons with opposite spins 
occupying the same site i. The effect of impurity is introduced in one of the lattice sites iA via repulsive 
interaction U a > 0. 

In the following we shall work in the Hartree-Fock approximation. To simplify the set equations we 
dropped the Hartree terms (Ui {n%a)), which means that we have only off-diagonal impurity induced disor- 
der in the system. So the corresponding Gorkov equation has the form: 



In the zero temperature the SC order parameter A$ and the total local charge m are given by self-consistent 
relations: 



1 f E i 



UiXi = -Ui- duJmG 1A (i,i\w + te), (3) 

7T ' 



— oo 



2 f Ef 



m = / dujlmG XL (i,i;u + ie), (4) 

J-oc 

where Ef denotes Fermi energy, and e is a small positive number. The local quasiparticle density of states 
(LDOS) can be written as follows 

AM =--lmG n (i > i;u + ie). (5) 

7T 

Eqs. l|2)-(|4} have been solved self-consistently using one site approximation at two dimensional lat- 
tice. The pairing amplitude has been found at the impurity site. The surrounding sites were treated 
homogeneously as in a bulk (clean) superconductor. 



3 Results and discussion 

It is well known that the local density of states at an impurity embedded in SC host has bound states 
within energy gap I^buik which are symmetrically located with respect to the Fermi energy. Especially this 
is true if the on-impurity Coulomb interaction has the same sign as SC one or if it is equal to zero. In this 
case bound states always have non-zero energies. However situation is quite different if the interaction has 
opposite sign to the SC one. One can show that under special conditions it is possible to get bound states 
which are exactly at the Fermi energy. 
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Fig. 1 Local density of states for Ubuik /W = —0.3 (a) and Ubuik /W = —0.6 (b). In Fig. la vertical and horizontal 
arrows show the change of LDOS with increasing interaction Ua- 



The example of the impurity local density of states (LDOS) is shown in the Fig.^where two situations 
are depicted. In Fig. ^ LDOS is plotted for Ubuik = —0.3 W, where W = 8t is a bandwidth, and for 
a number of on-impurity interactions Ua > 0. One can see that to get zero energy states it is not enough 
to have Ua of opposite sign to Ubuik- In this case bound states have non-zero energies for all Ua- If we 
increase Ubuik (see Fig. \T]p) and calculate LDOS one can find that those bound states do approach zero 
energy. It can be read from Fig. |2]where the positions of the bound states are plotted as a function of Ua 
for different values of the Ubuik- 

However it is very difficult to judge numerically if those states are really zero energy states. Therefore 
we supplemented our study by simple analytical calculations with constant bulk bare DOS. (Note that in 
numerical calculations we have chosen 2D tight binding DOS with Van Hove singularity in the middle of 
the band.) It turns out that the zero energy solution occurs for only one value of Ua which can be calculated 
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E/W 

Fig. 2 Position of the bound states as a function of Ua for Ubuik = —0.3 W (a) and Ubuik = —0.6 W (b). 
from 

A A ee Uaxa = -2nt 2 p(0), (6) 

where p(0) is the bulk bare density of states at the Fermi energy. 

To get such a large value of Aa we have to take large Ua as well as large value of Abuik because 
Xa in Eq. |6] depends on it via proximity effect. Physically we can imagine that there could be a strong 
on-impurity Coulomb repulsion. However as we can see in the Fig. |3]we cannot take large values of Ua 
because tends to a constant value as we increase Ua- The other possibility is to have large values of 




U A /W 

Fig. 3 Aa/ Abuik plotted as a function of Ua- 

Abuik which in fact is impossible to achieve in real materials. We need large values of the bulk pairing 
potential, much larger than the bandwidth. 
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4 Conclusions 

From the present calculations we conclude that it is in general not possible to get states at exactly zero 
energy in the system consisting of a single impurity embedded in BCS superconductor described by the 
negative U Hubbard model. However we were able to find one condition at which such situation can take 
place, which seems to be physically impossible to meet it in reality. Note the present calculations were 
performed assuming a rather small scattering rate e = 10 _3 < (Eq. [5] Fig. [2j. Of course taking a large 
enough e will be sufficient to broad and eventually connect the two neighbour bound states into one at zero 
energy. 

In this paper we have shown results for a half-filled band but we have also checked the situation for 
other fillings and as far as the bound states are concerned we got no qualitative differences. On the other 
hand present calculations are not fully self-consistent ones. We hope that proper treatment of the problem 
taking into account the Hartree term in Eq. |2] and calculating &buik at each surrounding site we will be 
able to give more precise conditions for 7r-like state appearance in such systems. 

Note also that relatively large values of the impurity potential U a (Fig. [H>) make Hartee-Fock approx- 
imation questionable so the results obtained here should be treated as qualitative. Our initial calculations 
show that fluctuations of the pairing parameter Aj around the central site would also influence the critical 
value of XJ a making it even higher. Such an analysis will be performed systematically for various band 
filings n. 
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